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ABSTRACT 

Context. Comparing models and data of pulsating stars is a powerful way to understand the stellar structure better. Moreover, such 
comparisons are necessary to make improvements to the physics of the stellar models, since they do not yet perfectly represent 
either the interior or especially the surface layers of stars. Because /3 Hydri is an evolved solar-type pulsator with mixed modes in its 
frequency spectrum, it is very interesting for asteroseismic studies. 

Aims. The goal of the present work is to search for a representative model of the solar-type star /3 Hydri, based on up-to-date non- 
seismic and seismic data. 

Mettiods. We present a revised list of frequencies for 33 modes, which we produced by analysing the power spectrum of the published 
observations again using a new weighting scheme that minimises the daily sidelobes. We ran several grids of evolutionary models 
with different input parameters and different physics, using the stellar evolutionary code ASTEC. For the models that are inside the 
observed error box of /? Hydri, we computed their frequencies with the pulsation code ADIPLS. We used two approaches to find the 
model that oscillates with the frequencies that are closest to the observed frequencies of /? Hydri: (i) we assume that the best model 
is the one that reproduces the star's interior based on the radial oscillation frequencies alone, to which we have applied the correction 
for the near-surface effects; (ii) we assume that the best model is the one that produces the lowest value of the chi-square (x^), i.e. 
that minimises the difference between the observed frequencies of all available modes and the model predictions, after all model 
frequencies are corrected for near-surface effects. 

Results. We show that after applying a correction for near-surface effects to the frequencies of the best models, we can reproduce the 
observed modes well, including those that have mixed mode character. The model that gives the lowest value of the^^ is a post-main- 
sequence model with a mass of 1.04 Mq and a metallicity slightly lower than that of the Sun. Our results underscore the importance 
of having individual frequencies to constrain the properties of the stellar model. 

Key words. Asteroseismology - Stars: solar-type - Stars: individual: /? Hydri 
1 . Introduction Il998h . its position in the Hertzprung-Russell (HR) diagram is, in 



ing the future of the Sun dPravins et al. 1993c b.a), making it a tion 6v to be about 55 //Hz, but we re unable to id entify individual 
particularly interesting object of study. mode frequencies. Subsequently, Bed ding et"aD ^07) observed 
Improvements to the fundamental parameters of y6 Hyi have Hyi during more than a week with the high-precision spec- 
been presented in a number of recent papers. Recent interfero- trographs HARPS and UCLES. Besides confirming the oscilla- 
metric measurements of /3 Hyi have yielded a n accurate (0.8%) tions detected in their previous observations in 2000, they were 
angulai- diameter for this stai- (North et al.l |2007) . Also, the able to identify 28 oscillation modes that included some mixed 
Hipparcos par allax of 13 Hyi has been i mproved f rom an uncer- modes of spherical degree / = 1 . Mixed modes occur in sta rs that 
tainty of 0.4% (iPerrvman & ESAll 1997 1) to 0.08% (Ivan LeeuwenI have l eft the main-sequence stage of their evolution (e.g. lOsakil 



l2oa 



ty or 0.4% (li^errvman & fctiAlliy^ / l) to U.USVc (Ivan l^eeuwenl nave l en me mani-seguence sialic oi Liieii evoiuuon i^e.g. iwsaKJi 
3). The combination of these two values gives a direct mea- 119751; lAizenman et al] |1977l) . and they provide useful informa- 



sure of p Hyi's radius with high accuracy. Moreov er, since the tion about the core. The presence of mixed modes, together with 
bolometric flux of this stai- is known (Blackwell & Lvnas-Gravl the availability of very precise non-seismic and seismic data for 

P Hyi, plac es the star in a priv ileged position for asteroseismic 



1. M. Brandao and G. Dogan contributed equally to this work. studies (e.g. lCuriha et al.ll2007h . 
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Fig. 1. Left panel: The position of (i Hyi in the HR diagram. The constraints on the fundamental parameters (Teff , L/Lq) are indicated 
by the 1-cr error box (solid) and on the radius by diagonal solid lines. We also show the corresponding 3-cr uncertainties by dashed 
lines. Two evolutionary tracks for the best models found by method 2 (cf. Table|5]) are plotted with dash-dotted and solid curves, 
representing the models with and without gravitational settling and diffusion, respectively. Right panel; the same as in the left panel 
but zoomed in. The selected models are marked by filled squares. 



Theoretical models of /3 Hyi based on its seismic and non- 
seismi c data have been published by Fernandes & Monteiro 
(2003), DiMauroetal. (2003), and Dog an et al.l (I2009D . 
[Fernandes & Monteiro (2003) examined the position of yS Hyi in 
the HR diagram by first considering the non-seismic data of the 
star In order to estimate the mass of /3 Hyi, they used avail- 
able seismic data, namely the large frequency separation, to re- 
move partially the helium-content vs mass degeneracy that exists 
when only non-seismic observational constraints are used. They 
also emphasized the usefulness of individual frequencies to con- 
strain the age of /3 Hyi due to the presence of mixed modes in its 
observed oscillation spectrum. Di Mauro et al. (2003) computed 
models of /5 Hyi, also based on its global parameters. They used 
the oscillation frequencies of /3 Hyi to compare with the model 
frequencies. Their theoretical models reproduced the observed 
oscillation spectrum of /? Hyiwell, as well as the observed large 
and small frequency separations, after they applied an ad-hoc 
shift to the computed frequencies. In fact, when comparing the 
computed and the observed frequencies, one should bear in mind 
that there may be an offset between them. This offset is well 
known from helioseismology and is also present when compar- 
ing the observed and computed frequencies for other stars. It 
arises from irnproper modelling of the surface layers of stars. 
iKjeldsen et al ] (l2008h used solar data to derive an empirical cor- 
rection for the near-surface offset which can be applied to stars 
other than the Sun. In our work, we apply this empirical correc- 
tion to the model frequencies of /3 Hyi bef ore comparing to th e 
observed ones. We extend the analysis of iDogan et al.l (l2009h . 
and also present a detailed discussion on the application of the 
near-surface correction. 



2. Observational constraints 

2.1. Non-seismic data 

T he most recent determ ination of the radius of /3 Hyi is given 
bv lKieldsen et al ] (12008ft . The radius was obtained by combining 



the in terferometric angu lar diameter of the star, 6 - 2.257+0.019 
mas dNorth et"al] l2007l ), with the revised H ipparcos parallax, 
TTp = 134.07 + 0.11 mas (Ivan Leeuwenll2007h . 

The luminosity, L, of a star can be obtained through the rela- 
tion. 



: 47rFboiC2/^p, 



(1) 



where Fboi is the bolometric flux and C is the conver- 
sion factor from parsecs to metres. To compute the lumi- 
nosity of /? Hyi we used the same value for the bolometric 
flux as North et all i2007f ). Fboi = (2.019 + 0.05) x lO"* W 
m"^ (jBlackwell & Lvnas-Grav 1998, the uncertainty on Fboi is 
from di Benedetto 19980 and the revised Hipparcos parallax 
(Ivan Leeuwen 200J) Adopting Lp = 3 .842 x lO^^W with an 
uncertainty of 0.4% dBahcall et al.ll200lh . we found L = 3.494 ± 
0.087 Lofor/JHyi. 

A number of determinations of the effective tempera- 



Favata et alJll997l:jBlackwell & Lvnas-Grav 


119981: di Benedetto! 


1998l:ISantos et al.l200.5l: Ida Silva et al.l2006 


Biazzo et al.l2007l: 


Bruntt et alj|2010l). We adopted the value of 


North etalJ (120071). 



which is derived from the direct measurement of the angular di- 
ameter. 

The mos t recent value for the metallicity of j3 Hyi is given by 
iBruntt etalJ ( l2010.) . [Fe/H] =^0^0 + 007. This value is in agree- 
ment wifli die one found bv lSantos et"al] (120 05), [Fe/ H1= - -0.08 + 
0.04 and with the one adopted bv lFernandes & Monteirol(l2003h . 
[Fe/ H1=-0. 12+0 . 047. In our analysis we adopted the metallicity 
from lBruntt et alj ( |2010|) . 

We calculated the mass fraction of metals, Z, from the metal- 
licity using the following approximation, valid for Population I 
stars which do not present the g-e lements enrichment seen in 
metal deficient stars dWheeler et al.l[T989l) : 



I. M. Brandao, G. Dogan et al.: Asteroseismic modelling of the solar-type subgiant star f} Hydri 



3 



Table 1. Stellar properties of p Hyi. The luminosity, L, and ra- 
dius, R, are expressed in solar units. stands for the angular di- 
ameter, TTp for the Hipparcos parallax, Teff for the effective tem- 
perature, [Fe/H] is the metallicity, and ZjX is the mass ratio of 
heavy elements to hydrogen. 



mode frequencies with coiTect I identifica tions is somewhat sub - 
jective. We followed the same approach as ^Bedding et alJ(l2010h . 
which involved using the ridge centroids as a guide and averag- 
ing multiple peaks into a single value. The remaining unidenti- 
fied peaks in the power spectrum are listed in Table[3] 





Value 




Reference 


6 (mas) 


2.257 ± 


: 0.019 


North et al. (2007) 


Tip (mas) 


134.07 


±0.11 


van Leeuwen (2007) 


^/Rg 


1.809 4 


: 0.015 


Kieldsen et al. (2008) 




3.494 4 


: 0.087 


Current work 


r.ff (K) 


5872 ± 


44 


North et al. (2007) 


[Fe/H] 


-0.10 H 


fc0.07 


Brunttetal. (2010) 


zix 


0.019 4 


0.003 


Current work 



,Pem..log(f)_.,og(|)_-l„.(f)^-l„g(|)^ 

where [Fe/H]s is the star's metallicity; Zpe and X are the iron and 
hydrogen mass fractions, respectively; and (Z / X)pi is the ratio for 
the s olar mixture. We used (Z/X)o = 0.0245 (iGrevess e & Noeli 

[1993 ). This gives (ZIX) = .019 ± 0.003 for ^6 Hy i 

From spectral analysis, iDravins & Nordlundl (11990') found 
vsin; = 2 ± 1 km s ' for fi Hyi. More rece ntly, [ Bruntt et ak 
('2010') found vsin/ = 2.7 + 0.6 km s ', and "Hekker & AertsI 
(^10) from spectroscopic line-profile an alysis, fou n d v sin ; = 
4.3 km s From their analysis, Hekk er & AertsI (1201 Oh at- 
tempted to determine the inclination angle, of fi Hyi, suggest- 
ing a value of 55 + 17° for this star. Thus the effect of rotation 
on the modelling of the structure of the star can be neglected. 
Similarly, since the resulting rotational splitting is comparable 
with the error in the observed frequencies (see below), in the 
present analysis we neglect the effects of rotation on the fre- 
quencies. 

The position of (3 Hyi in the HR diagram is shown in Fig.[T] 
and the fundamental parameters we adopted are given in Table[T] 



Table 2. Observed oscillation frequencies in p Hyi (in //Hz) re- 
sulting from the revised analysis, listed in ascending radial order 
within each column. The rows are in ascending /, and each row 
includes frequencies within Av-sized-bits of the frequency spec- 
trum, is used for the modes whose S/N was too low for a 
clear extraction. 



/ = 


/= 1 


1 = 1 


/ = 3 


660.74 ± 2.43 








716.68 ±3.00 




711.24±2.I3 




774.79 ± 2.20 


802.74 ± 1.69 


769.97 ± 0.99 


791.66 ± 1.35 


831.86 ±2.43 


857.32 ±0.86 


825.86 ± 1. 1 8 




889.15 ± 1.23 


912.91 ±0.86 


883.35 ± 0.89 




946.11 ±0.91 


959.98 ±0.89 
987.08 ± 0.87 


939.97 ± 0.97 




1004.32 ±0.86 


1032.99 ± 0.86 


999.40 ±0.91 




1061.66 ±0.95 


1089.87 ± 0.88 


1057.00 ± 0.86 




11 18.67 ±0.88 


1147.35 ±0.91 


1115.20 ± 1.06 




1177.76 ±0.97 


1203.54 ± 1.01 


1172.98 ±0.86 


1198.16 ± 1.23 


1235.31 ± 1.09 


1320.42 ±0.94 
1378.92 ± 1.39 







Table 3. Unidentified observed peaks with S/N > 3.5. 



V (/iHz) 



753.12 + 1.57 


1013.42 ± 


1.50 


1130.36 ± 


1.30 


828.70 ± 1.83 


1025.80 ± 


1.68 


1134.32 ± 


1.63 


845.02 ± 1.61 


1037.90 ± 


1.63 


1167.62 ± 


1.10 


868.60+ 1.13 


1065.12 ± 


1.59 


1256.78 ± 


1.60 


911.88 ± 1.76 


1070.00 ± 


1.43 


1383.20 + 


1.75 


1010.20 ± 1.91 


1084.20 ± 


1.57 


1462.62 ± 


1.92 



2.2. Seismic data 

Asteroseismic observ ations of /3 Hyi have been reported by 
iBedding et al.l (l2007h . They found an excess power centred 
around 1 mHz with a peak amplitude of ~50 cm s and the os- 
cillation frequencies show a comb-like structure typical of solar- 
like oscillations with a large frequency separation of the / = 
modes, Avo, of 57.24 ± 0.16 //Hz. They also identified 28 mode 
frequencies in the range 0.7 < v < 1 .4 mHz with the angular de- 
grees I -0, 1, and 2, three of which were identified as / = 1 mixed 
modes. In this work, we used the updated list of 33 observed fre- 
quencies given in Table|2l T o produce these, we reanalysed the 
2005 dual-site observations (iBedding et al.ll2007h using revised 
weights that were adjusted using a new method that minimises 
the sidelobes ( H. K jeldsen et al., in preparation). This method is 
described by iBedding et al.l (l2010l), who applied i t to multi-site 
observations of Procyon (see also lArentoft et al] 12009) . In the 
same way as for Procyon, we extracted oscillation frequencies 
from the time series of /3 Hyi using the standard procedure of 
iterative sine-wave fitting. The finite mode lifetime causes many 
modes to be split into two or more peaks which, coupled with the 
presence of mode bumping, meant that deciding on a final list of 



3. Modelling Hyi 

3.1. Input pliysics to tlie models 

To compute the evolutionary models we used the Aarhus STellar 
Evolution Code', ASTEC (Christensen-Dalsgaard 2008a). The 
following assumptions were made: spherical symmetry, no ro- 
tation, no magnetic field and no mass loss. In the computa- 
tion we used the OPAL 2005 equation o f state tables (see 
Rogers & Navfonovl l2002h . OPAL opacities ("Iglesias & RogersI 
1996) complemented by low-temperature opacities from 
Fer ^son et al.l (120051) . the solar mixture from lGreves se & Noels 
(1I993I) . and the NACRE nuclear reaction rates dAngulo et al. 
Il999h . We considered an atmospheric temperature versus op- 
tical depth relation which is a fit to the quiet-sun relation of 
Vernazza et al. ( 1976). 

Convection was treated according to the stand ard mixing- 
length theory (MET) from lBohm-Vitens3 (Il958l) . where the 
characteristic length of turbulence, /mix, scales directly with the 
local pressure scale height. Hp, as l^i^ - aHp, leaving the scal- 
ing factor a as a free parameter. MET makes use of another 
free paramete r associated with the am ou nt of the core overshoot- 
ing, Q-Qv. Both lDi Mauro et an(l2003l) and lFernandes & Monteird 
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(l2003h found that models at the position of p Hyi in the HR dia- 
gram are not affected by convective overshooting, so we decided, 
for this work, not to consider it in our models. 

Diffusion and s ettling were treate d in t he approxima- 
tions proposed b y iMichaud & Proffitt" ("19931). We refer to 
IChristensen-Dalsg aard (2008a) for a detailed explanation. 

3.2. Grids of models 

We calculated two grids of evolutionary tracks. Grids I and II, 
with the input parameters shown in Table|4] In Grid II we in- 
cluded diffusion and gravitational settling of helium. 

For each grid, we took those models whose parameters 
were within the 3-cr uncertainties derived from the obser- 
vations of /3 Hyi, and computed the corresponding oscilla- 
tion frequ encies using the Aarh us adiabatic oscillation code 
ADIPLS (IChristensen-DalsgaardI r2008bl) . The theoretical fre- 
quencies were then compared with the observed ones in order 
to find the best-matching model. 

4. Near-surface corrections 

It is well known from helioseismology that there is a systematic 
offset between the observed and the computed oscillation fre- 
quencies of the Sun. This offset, which is nearly independent of 
the angular degree , I, of the mode and affects the highest frequen- 
cies the most (Chr istensen-Dalsgaard & Thomp son 1997), arises 
from the improper modelling of the surface layers. Therefore, 
the offset is also expected to be present when comparing the ob- 
serve d and computed frequencies for other stars. Kields en et al.l 
(l2008i) used the solar data to derive an empirical correction for 
the near-surface offset, which can, in principle, be applied to 
other stars. 

For the Sun, the difference between the observed, Vobs, and 
computed frequencie s of the best model, Vbest, was shown by 
iKieldsen et al.l (|2008|) to be well approximated by a power law 
fit given as 



Vobs(n) - Vbest(n) = a 



Vobs(n) 



Vo 



(3) 



where n is the radial order, and vo is a reference frequency that is 
chosen to be the frequency at maximum power Since the offset 
is independent of I, the authors considered only radial (/ = 0) 
modes. Note that the 'best model' is the one that best represents 
the interior but still fails to model the near-surface layers. They 
also argued that the frequencies of a reference model, Vref , which 
is close to the best one, can be scaled from Vbest by a factor r, i.e.. 



Vbest(«) = '"Vref(n)- 



Then Eq. (O becomes 



Vobs(«) - '"Vref(n) = a 



Vobs(n) 



Vo 



(4) 



(5) 



IKieldsen et all (l2008l) . using the data and models for the Sun, 
found b = 4.90. Using this value for b it is possible to deter- 
mine r and a from Eqs (6) and (10) in Kieldsen et al. (2008)). 
Assuming that a similar offset occurs for other solar-like stars, 
they showed how to use the solar b to determine r and a from 
the observed frequencies of I - Q modes, and, consequently, use 
Eq. (|5]l to calculate the correction that must be applied to the fre- 
quencies computed for all I for a given stellar model. They noted 
that the correction applied to the mixed modes should be less 
than that for the pure p modes (see the next section). 



5. Selecting the best model 

We considered two methods to find the model that to the closest 
extent possible reproduces the observed non-seismic and seismic 
data for/? Hyi. For all models whose parameters were within the 
3-cr uncertainties derived from the non-seismic observations of 
P Hyi, we calculated the r and a values, following lKieldsen et al.l 
(l2008h . using b = 4.90 and vo = 1000 yuHz. 

In the firs t meth od we followed closely the work of 
IKieldsen et al.l (|2008|) and considered the best representative 
model to be the one having the value of r closest to 1, which 
means, from Eq. 



Vref(n) ~ Vbest(n)- 



(6) 



Using the values of r and a found for this model, we then com- 
puted the correction factor to be applied to the model frequencies 
and compared them with those observed. Note, however, that in 
this method r is calculated using only the observed radial modes 
together with the corresponding model radial modes. If we as- 
sume that the best model is the one which has r closest to 1 
then we are assuming the model that has the radial frequencies 
matching most closely the observed / - modes also has the 
nonradial frequencies matching best the observed I - 1 and 2 
modes, which may not be true. 

In the second method we performed a statistical analysis, us- 
ing the test. Our goal was to find the reference model that, 
after applying the surface corrections, had the individual fre- 
quencies closest to the observed ones, for all I, i.e. we found 
the minimum of 



X 



Vi-ef,corr(n, ' Vobs(n, 



0-(Vobs(«, 0) 



(7) 



where is the total number of modes considered, Vref,corr(n, 
are the frequencies of modes with radial order n and degree I 
of a reference model, corrected for the surface effects, and cr 
represents the uncertainty in the observed frequencies. 

The correction term, as shown in the right-hand side of 
Eq. ([3]) can only be applied to the frequencies of the best model. 
In order to compute the correction term, since we have a set 
of reference models and we do not know which one is the best 
model, we assume that the corrected best Vbest.coiT and reference 
model Vi-ef,con- frequencies also scale as 



Vbest.con" — ^l^ief,coiT (n,0). 



(8) 



We note that this is a good approximation because the surface 
corrections to the frequencies are much smaller than the frequen- 
cies themselves. Moreover, since Vbest.con- - Vobs we have 



Vref,corr(«, 0) = Vref («, 0) + j 

Thus, Eq. (|7]l becomes 



Vobs(«, 0) 



Vo 



(9) 



X 



nJ 



'y,ef(«,/) + (f)[^f -Vobs(«,0 



o"(Vobs(n,0) 



(10) 



Note that, in practice, the term (fl/r)[vobs(n, 0)/vo] on the right- 
hand side of Eq. (|9]l was replaced by (fl/r)[vi-ef(n, 0/vo]*- The 
reason is to enable us to correct all the reference model frequen- 
cies, instead of only the frequencies having the same radial order 
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Table 4. Parameters used to compute the evolutionary tracks. M/Mq is the mass in solar units, Z/X is the initial ratio of heavy 
elements to hydrogen abundances, and Y the helium abundance. 



Parameter 


Grid I 




Grid 


11 




1.00- 


1.18 (with steps of 0.02) 


1.00 


- 1.18 (with steps of 0.02) 


Z/X 


0.010 


- 0.030 (with steps of 0.004) 


0.010 - 0.030 (with steps of 0.004) 


Y 


0.24- 


0.30 (with steps of 0.02) 


0.24 


- 0.30 (with steps of 0.02) 


Mixing length 










parameter (a) 


1.4-2 


.0 (with steps of 0.2) 


1.4- 


2.0 (with steps of 0.2) 


Diffusion & 










gravitational settling 


None 




He 





as those observed. The difference between these terms is mini- 
mal for the best-fitting model, and so we are safe in performing 
the replacement. 

To compute we used all the observed and computed 
/ = 0, 1, and 2 frequencies. However, as already mentioned, 
some of the observed 1=1 modes are mixed modes and the scal- 
ing is not valid for them. Mixed modes should not be affected by 
the surface layers as much as the p modes (Kieldsen et al. 2008), 
so the coiTection term should be small for the mixed modes. 
Specifically, at a given frequency we expect the near-surface ef- 
fects to scale inversely with the mode inertia, which is much 
higher for the mixed modes than for the p modes; thus we scaled 
the correction term by the inverse of the ratio Q„i between the 
inertia of t he mode and the i nertia of a radial mode of the same 
frequency dAerts et alj|2010() . Taking that into account, Eq. (fTOl l 
becomes 



X 



-y 

n,l 



fv..eK«,/) + (^)(f)[^f-Vo..(«,0 



0-(Vobs(«, 0) 



(11) 



In this method, our best model is the one that minimises Eq. ( fTTT i. 

Contrary to the first method, the test takes into account 
all the individual frequencies, including mixed modes. This test 
also includes the uncertainties in the observations. Thus, it is 
possible that the model with the lowest does not have the r 
that is closest to 1 . 



6. Results and Discussion 

The parameters of the best models found for Grids I and II, us- 
ing the two methods described in Sec. |5j are shown in Table|5] 
Figures|2]and[3]show the echelle diagrams for/? Hyi. An echelle 
diagram shows the frequency spectrum divided into segments 
equally spaced by the large frequency separ ation, afte r these seg- 
ments are stacked in the vertical direction jCirec et al. 1983) . In 
these figures the observed frequencies of yS Hyi are compared 
with the theoretical frequencies of the best models from Grid I 
(upper panel) and from Grid II (lower panel), both before (left 
plot) and after (right plot), applying the near-surface corrections. 
The model frequencies are represented by open symbols and the 
observed frequencies (cf. Table|2]l by solid symbols, while the 
asterisks represent the unidentified peaks (cf. TableO, which 
may correspond to genuine modes, sidelobes, or noise peaks. 
The relative sizes of the open symbols reflect the expected mode 
amplitudes dChristensen-Dalsgaard et all Il995t) . The so-called 
mixed modes reveal themselves in the echelle diagrams, break- 
ing the regularity of the ridges. The models predict mixed modes 
with all nonradial degrees, however mostly with too small ampli- 
tudes to be observed. On the other hand, some of the observed 
modes match well the mixed modes with / = 1 (see, e.g., the 



right panels of Fig. [3] If we inspect Figs.|2]and[3] it is clear that 
the agreement between the observed and model frequencies is 
much better when method 2 is used (Fig.O. This is due to the 
fact that in this method all the available seismic constraints were 
involved in selecting the best model. 

It is seen from Table|5]that the model that has the r value 
closest to unity does not produce the lowest value. The model 
with the lowest still has an r satisfactorily close to unity. So, in 
addition to finding a model that represents the stellar interior rea- 
sonably, method 2 makes sure that all the available seismic con- 
straints are simultaneously reproduced and so the fit, and hence 
the accuracy of the model, is improved substantially. This shows 
the importance of the individual modes in constraining the range 
of models to represent the observed star. Mixed modes, in par- 
ticular, put strong constraints on the model properties, especially 
on the evolutionary stage. For instance, we can see from the right 
panels of Fig. |2] that two models resulting from method 1 have 
the two highest values due to failing to match particularly the 
observed I = 1 mixed modes. In the upper right panel, the model 
is too massive and it matches the rest of the seismic constraints 
before it is sufficiently evolved to have mixed modes, whereas 
the model in the lower right panel does have mixed modes, al- 
though the predicted mixed modes do not match those observed. 

In gener al, we found that the empirical surface corrections 
proposed by 'Kieldsen et al.l (l2008h work very well for yS Hyi as 
seen from Figs.|2]and|3] although there is room for improvement 
for high-frequency modes of / = 1 . The reason for the subop- 
timal agreement for those modes is that the correction term is 
determined using only the / = observed modes, whose fre- 
quencies span a smaller range than those of the Z = 1 modes. 
Thus, radial modes with higher frequencies need to be detected 
in order to improve the agreement for the higher frequency / = 1 
modes. Note that the change in the large frequency separation of 
the models after applying the near-surface coiTection is around 
0.8 /vHz, which is larger than the given uncertainty of the ob- 
served large separation. This should be taken into account when 
modelling through a pipeline analysis that uses the large sep- 
aration as input. It is encouraging to see that we can observe 
/ = 3 modes, and that some of the unidentified modes are also 
close to the model frequencies, namely 753.1 yuHz (/ = 1?) and 
1462.6 yuHz (/ = 0?). 

Even though the best model seems to be the one without dif- 
fusion, we do expect that within a star diffusion occurs. The two 
selected models, the ones with and without diffusion, resulting 
from the method 2 are in fact compatible and both could be fur- 
ther fine-tuned. 



7. Summary and Conclusions 

We computed two grids of evolutionary models, using the 
ASTEC code, in order to find the model that most closely repro- 
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Table 5. The parameters of the best models found for Grid I (no diffusion) and II (He settling and diffusion), for each of the two 
methods. See text for details on the methods. The mass, M, luminosity, L, and radius, R, are expressed in solar units. T^a is the 
effective temperature, Y and Z are the initial helium and heavy-element abundances, [Fe/H] is the metallicity at the surface, and a 
is the mixing-length parameter Also r and a are factors used to compute the correction term, Avq/, and Avoa are, respectively, the 
large frequency separation before and after applying the surface correction to the model I - modes. The values of;^'^ are those 
calculated after correcting the frequencies for the near-surface effects. 





Grid! 


Grid II 


Pin rn m ptpr 


Method 1 Method 2 


iVlcLIlUU i 


iVlCLIlUU Z 


M/Mo 


1.16 


1.04 


1.04 


1.04 




1.832 


1.785 


1.790 


1.786 


LIU 


3.433 


3.485 


3.432 


3.338 


r.ff (K) 


5810 


5908 


5877 


5843 


Age (Gyr) 


4.705 


6.114 


7.390 


7.297 


Z 


0.0204 


0.0124 


0.0075 


0.0075 


Y 


0.30 


0.30 


0.24 


0.24 


[Fe/H] 


0.088 


-0.133 


-0.416 




a 


1.4 


1.8 


2.0 


1.8 


r 


1.0000 


0.9995 


1.0000 


1.0009 


a (juHz) 


-4.80 


-3.14 


-2.43 


-3.11 


Avoi ifMz) 


58.977 


58.488 


58.243 


58.400 




57.678 


57.652 


57.600 


57.577 




19.086 


1.183 


26.226 


2.642 



duces the seismic and non-seismic data of ji Hyi. The parameters 
used for each grid are given in Table[T] We computed the oscilla- 
tion frequencies for the models that lie inside the 3-cr error box 
of the star's position in the HR diagram with the ADIPLS code, 
and compared them with the observed frequencies. There is an 
offset needed to be taken into account in this comparison, due to 
improper modelling of the near-surface layers of the stars. We 
used the approach proposed by Kjeldsen et al. (2008) to correct 
the computed frequencies of fi Hyi from this offset. We used two 
methods in order to find the model that reproduces the observed 
oscillation frequencies of fi Hyi. 

In our analysis, we argue that the method involving the 
test, method 2, is the most robust way to find the best model, 
since it takes into account all the individual frequencies, the 
mixed modes, and also the uncertainties on the observed fre- 
quencies. Analysing the echelle diagrams of the representative 
models found with method 2 (cf. Sect.|6l), we see that the surface 
correction works very well for / = modes, and for I - 1 and 
2 modes with frequencies lying in the frequency range of the 
observed radial modes. This was expected since the correction 
term was computed using only those radial modes. Observed 
/ = modes with higher frequencies are thus needed in order 
to improve the surface correction. 

Our best models give M=1.04 and an age of 6.1 - 7.3 
Gyr for fi Hyi, depending on the inclusion of gravitational set- 
tling and diffusion of helium. In either case, the radius is found 
to be 7? ~ 1 .785 Rq, which is in good agreement with the one 
determined by interferometry, R - 1.809 + 0.015 Rq. However, 
there are other models fitting the data similarly well. We used 
the parameters of those models (with;^f^ < 10) to determine the 
internal error regarding our analysis. We calculated the mean 
value, and the uncertainties were taken as the standard devia- 
tion. We found M = 1.08 + 0.03 Mq, age = 6.40 + 0.56 Gyr, 
and R - 1.811 ± 0.020 Rp. These resul t s are also consistent 
with the results of iFernandes & Monteird (12003 ^. who derived, 
M = 1.10!° °^ Mq and M = 1.09 ± 0.22 Mq, through the HR di- 
agram analysis and Avq, respectively, and a stellar age between 
6.4 and 7.1 Gyr. 
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Fig. 2. Echelle diagrams for /? Hyi, with a frequency separation of Av - 57.5 yuHz, before (left plot) and after (right plot) application 
of the near-surface corrections to the model frequencies. Shown are the frequencies of the selected models using method 1, when 
including no diffusion (upper panel) and diffusion (lower panel). In method 1, the best model was selected using the radial {I - 0) 
modes alone (see the text for details). The solid symbols show observed frequencies (Table|2]i, asterisks the unidentified peaks 
(TableO, and the open symbols the model frequencies. Circles are used for I - modes, triangles for Z = 1, squares for I - 2 and 
diamonds for I -?>. Open symbols are scaled to represent the relative amplitudes of the modes as predicted by the models. 
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Fig. 3. The same as Fig.|2]but for the best models without (upper panel) and with (lower panel) diffusion, selected using method 2, 
which takes into account the observed and identified modes with all degrees available. 



